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I make a number of observations concerning discretization errors in twisted-mass lattice QCD that 
can be deduced by applying chiral perturbation theory including lattice artifacts. (1) The line along 
which the PCAC quark mass vanishes in the untwisted mass-twisted mass plane makes an angle 
to the twisted mass axis which is a direct measure of 0(a) terms in the chiral Lagrangian, and is 
found numerically to be large; (2) Numerical results for pionic quantities in the mass plane show the 
qualitative properties predicted by chiral perturbation theory, in particular an asymmetry in slopes 
between positive and negative untwisted quark masses; (3) By extending the description of the "Aoki 
regime" (where m q ~ a 2 Aq CD ) to next-to-leading order in chiral perturbation theory I show how 
the phase transition lines and lines of maximal twist (using different definitions) extend into this 
region, and give predictions for the functional form of pionic quantities; (4) I argue that the recent 
claim that lattice artifacts at maximal twist have apparent infrared singularities in the chiral limit 
results from expanding about the incorrect vacuum state. Shifting to the correct vacuum (as can be 
done using chiral perturbation theory) the apparent singularities are summed into non-singular, and 
furthermore predicted, forms. I further argue that there is no breakdown in the Symanzik expansion 
in powers of lattice spacing, and no barrier to simulating at maximal twist in the Aoki regime. 



I. INTRODUCTION 

Twisted mass fermions p| at maximal twist H H are presently under intensive study both with simulations and 
theoretical methods. 1 The standard application involves two degenerate quarks, and this is the theory considered 
here. 2 The long distance physics (vacuum and pionic properties) of twisted mass lattice QCD with two degenerate 
flavors (tmLQCD) can be studied using the methods of effective field theory, allowing a systematic dual expansion 
in the quark mass, m q , and the lattice spacing, a, in which non-perturbative effects are parameterized by a priori 
unknown low energy constants (LECs) d U EH . The resulting theory is called twisted-mass chiral perturbation 
theory (tmxPT) [Tl|. Its consequences in the vacuum and pionic sector have been worked out to next-to-leading 
order (NLO) in the regime with m q > aA 2 (the "generic small mass" or GSM regime 3 ) [l2L IT^| , and at leading order 
(LO) in the regime with m q — a 2 A 3 (the "Aoki" regime) El El El El EH < Here A is a scale of the order of the 
QCD scale, A QC d- 

The present paper is a collection of observations that follow from tnr\;PT concerning discretization errors in tm- 
LQCD. On the practical side, I point out in seclTTIsome new applications of results contained in Ref. EH ■ In particular, 
I describe how the observed difference between two methods of determining maximal twist that are being used in 
simulations can be understood using tm%PT. This difference provides a direct measure of the size of discretization 
errors, or equivalently a measure of the LECs related to such errors. I also point out the relationship of these methods 
to a third (theoretically interesting but apparently less practical) method proposed in Ref. Ell- 
in sec. IIIII I note that the dependence of physical quantities on the untwisted quark mass at fixed twisted quark 
mass gives another set of measures of the size of discretization errors, some of which are related by tm%PT at NLO. 
These results are contained in Ref. [lIj but have not been compared to the available numerical data, and I think it 
is useful to provide illustrative plots of the expected types of dependence. 

When making these plots the question arises as to how to extend them into the Aoki regime, which is the regime in 
which there are phase transitions caused by competition between quark mass effects and those due to discretization 
errors. Previous studies, in Refs. El El El El El, 

worked to LO in this regime, and here I extend these results to 
NLO. This extension, described in sec. II VI turns out to require only one non-trivial additional low-energy coefficient. 
The results allow one to see how the phase transition lines, which are linear at LO, can become curved at NLO, 
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1 For recent reviews see Refs. 

2 Extension to non-degenerate quarks has been explained in Ref. 0. 

3 Note that the GSM regime is here defined to include the region in which m q ^> aA 2 , since the NLO expressions remain valid in this 
region, although the terms proportional to a 2 will become of NNLO. Of course, for chiral perturbation theory to apply one must always 
have that m q <C A. In practice, for present lattice spacings, this latter condition implies that m q does not exceed aA 2 . 
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and similarly how the lines of maximal twist using different definitions behave. These results may be amenable to 
numerical investigation, but, in any case, are of theoretical interest in understanding the relationship between different 
definitions of twist angle and critical mass. 1 also point out that, while a classical analysis is sufficient when working 
to NLO in the Aoki regime, extension to next-to-next-to-leading order (NNLO) requires the inclusion of loop effects. 

My final topic, discussed in sec.E] is the conclusions of Ref. |17|. ft is argued in [l7| that, if one makes an estimate 
of the critical mass that has an error of size aA 2 , then, when working at maximal twist, physical correlation functions 
contain apparently infrared (IR) divergent discretization errors proportional to (a / 'TO 2 ) 2fe , with k > 1 an integer. 
Such errors, while consistent with automatic O(a) improvement, would restrict one to working with m q > aA 2 , and 
suggest a breakdown of the Symanzik expansion 18] for smaller quark masses. It is then argued that, either by using 
an "optimal" choice of critical mass, or using a non-perturbatively improved quark action, these IR divergent errors 
at maximal twist can be reduced in size to a 2 (a 2 /m 2 . ) fe . In this way the bound on the quark mass is reduced to 



> a 2 A 3 , i.e. one can work in the GSM regime but not in the Aoki regime. 



My observation is that the IR divergences are an indication of expanding around the wrong vacuum in the presence 
of large perturbation. Using tm^PT one can sum these corrections to all orders and obtain a finite and smooth 
prediction for the dependence of physical quantities on the underlying untwisted and twisted quark masses. Indeed, 
this is precisely what was done previously in Refs. |l2lll3| |. I argue further that there is no breakdown of the Symanzik 
expansion, suitably interpreted, and no barrier to working at maximal twist in the Aoki regime. 

These considerations do not, however, invalidate the conclusion of Ref. that one should use an appropriate 
"optimal" choice of critical mass to reduce the size of discretization errors. This was proposed previously using 
tmxPT in Refs. 0,E1> an d the criteria for determining the optimal choice agree. My observation here is that, since 
tmxPT is built upon the Symanzik expansion used by Ref. |17| . the result concerning optimal critical masses had 
to agree. In fact, tm^PT goes beyond the Symanzik expansion by adding non-perturbative information concerning 
spontaneous chiral symmetry brea king . The only additional assumption that one has to make is that the effective 
field theory framework is valid \19L |2£| . 

It is useful to keep in mind during the subsequent discussions the physical values of quark masses that correspond to 
the different regimes. If one takes a -1 = 2 GeV (a relatively small lattice spacing for present dynamical simulations) 
and A « Aqcd ~ 300 MeV, then aA 2 (the size of quark masses in the GSM regime) and a 2 A 3 (the size in the Aoki 
regime) are about 45 and 7 MeV respectively. Since we are interested in calculating with light quark masses varying 
from m s trangc ~ 80 — 100 MeV down toward mij g ht ~ 3 — 4 MeV, most simulations will be done in the GSM regime, 
with a significant tail entering the Aoki regime. Thus it is important to understand both regimes in detail. 

In this paper I work in the "twisted basis" in which the tmLQCD action is 



+ V M )-^V*V A1 + m + z 7 5T3Mo]^^) I ■ (1) 



Here tpi is the dimensionless bare lattice field, fj,o is the bare twisted mass and too the bare untwisted mass. The 
coefficient of the Wilson term, r, is typically set to unity in simulations; the considerations that follow hold, however, 
for any O(l) choice of this coefficient. Maximal twist is obtained by setting too to an estimate of its critical value 
to c (which, as will be seen, can depend on /j,q). The determination of to c is a key issue in tmLQCD, and the topic of 
much of the discussion which follows. I only note here that this is a non-trivial issue as no symmetry is restored when 
too = m c . By contrast, setting fiQ = restores parity and flavor symmetries, and physical quantities are symmetrical 
under /xo — > — Mo- 

I use the twisted basis, rather than the more continuum-like "physical basis", used for example in Ref. |17|. since 
the twisted basis is that usually used in simulations, allowing a more direct connection to the numerical data. It is 
also the basis used in the tm^PT calculations of Ref. [lj] , which I use extensively. 

The physical twisted and untwisted masses are [i = Z^^o/a and to = Z m (mo/a — m c ), respectively, with = 1/Zp 
and Z m — 1/Zg renormalization constants. 4 The physical quark mass is then m q = \J /i 2 + m 2 , and to work at 
maximal twist requires setting to = (i.e. setting too = am c ). The accuracy with which to c must be determined in 
order to maintain automatic 0{a) improvement depends on how close to the chiral limit one is working. If m q 3> aA 2 , 
then discretization errors are small, and the error in m c can be of 0(a). 5 This possibility is part of the topic of sec.lVl 
In the GSM regime, one must incorporate O(a) discretization effects into m c 0,E1> an d thus reduce the error to 
0(a 2 ). To work in the Aoki regime and maintain automatic 0{a) improvement requires the uncertainty in m c to be 



4 Note that I define m c to be dimensionful. 

5 Here, and in the following, I often use a notation in which appropriate factors of A are implicit. 
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further reduced, to 0(a 3 ) [H[l3|. Following Ref. , I refer to choices of critical mass which are accurate to 0(a 3 ) as 
"optimal", and denote them by m' c (and the resulting untwisted quark mass by ml). I discuss these choices in secs.HTI 
and IIIII One noteworthy choice of m Cl which, as stressed in Ref. [l^ . is not optimal, is given by the position where the 
pion mass vanishes (an end-point of the Aoki phase, if present). This differs from optimal choices by 0(a 2 ). Finally, 
if one works at NLO in the Aoki regime, it is natural, as shown in sec. IIVI to incorporate an 0(a 3 ) shift into m c . I 
denote the resulting critical mass and untwisted physical quark mass by m" c and m", respectively. The uncertainty 
in m" turns out to be of 0(a 5 ). The fact that only odd powers of a appear has been noted also in Ref. [l7j . 

So as not to clutter the notation with too many "primes" , I use the same (unprimed) symbol for m q in each regime, 
although it is defined using the untwisted quark mass appropriate to the given regime. Thus, for example, in the Aoki 
regime, m q = -\/m 2 + m" 2 . 

The variables mo /a and fio/a (or m and fi, or m! and fj,, etc.) map out what I will refer to as the "mass plane". 
The origin in this plane is defined to be where /io = and m /a is set to the critical value appropriate to the regime 
of interest. With the definitions given above, m q is then the distance from the origin. Once m c is estimated, one 
possible definition of the twist angle (that which relates most directly to that used in the continuum) is the angle 
relative to the positive untwisted mass axis. This angle depends on the choice of critical mass, but I refer to it as lvq 
in all regimes. It is most useful in the GSM regime, where it is defined by 

uj Q = tan (/i/m') . (2) 

In the Aoki regime, at NLO, one replaces m! with m" . 

Finally, I collect the notation I use for low energy coefficients (LECs) in tm^PT, largely taken over from Ref. 0] . At 
leading order, there are two continuum LECs: /, the decay constant in the chiral limit normalized to = 132 MeV, 
and Bo, defined so that m 2 — 2Bom q . For quark masses renormalized at 2 GeV in the MS scheme, Bq w 2.5 GeV. It 
is then convenient to introduce ft = 2Bo/i and rh = 2B m and M = jl 2 + m 2 so that, at LO, m 2 = M. At NLO 
in the continuum, there are the standard Gasser-Leutwyler coefficients, the Li [5lJ. Discretization errors introduce 
other LECs. At LO in the GSM regime, there is a single constant Wq, which determines the 0{a) shift in m c Hi- In 
particular, defining a — 2W^a, then the shifted untwisted mass is rh' = 2Bom' = rh+a, and = M' = \/ fi 2 + (m') 2 . 

Discretization errors of NLO in the GSM regime are determined by the dimensionless constants W, W\o, W, and 
It is convenient to introduce the following quantities: 



. _ 16a(W + Wio/4) . 16a(W + W w /2) , _ 16a 2 W 

6 W = j 2 , d w = j 2 , w = ^ 2 . (>») 



Syy and give dimensionless measures of the size of discretization errors in physical quantities, while the mass- 
squared w' appears additively in the result for the pion mass-squared. Note 5y/ and 8^, are linear in a, while w' is 
quadratic. Further LECs enter at NLO in the Aoki regime, and will be defined in sec. IIVI 

II. COMPARING DEFINITIONS OF MAXIMAL TWIST 

As noted in the introduction, one must use a non-perturbative method of sufficient accuracy in order to determine 
when one is at maximum twist. This section concerns the relation between two methods being used in present 
simulations [which I refer to as methods (i) and (ii)] and their relationship to a third method [method (iii)] proposed 
in Ref. 01 ■ Throughout this section I assume that the quark masses are in the GSM regime. 

My observations are summarized in Fig. The line in the mass-plane defined by method (i) (in which the PCAC 
mass is set to zero) makes an angle 5u> to the twisted mass direction [and thus, by construction, to the line defined 
by method (ii)]. This angle is proportional to the lattice spacing, and to one combination of the LECs describing 
discretization errors. It can be determined from present simulations, and the results show that this discretization 
error is not as small as one might hope, since it is determined by a scale A w 0.7 GeV rather than by A sa Aqcd- 
I also observe, following Ref. |13|. that method (iii), previously argued to be difficult to implement in practice, is 
actually equivalent to method (ii) in the GSM regime (and thus has already been used). 

Both method (i) and (ii) are based on the continuum result that, at maximal twist, charged axial currents in the 
twisted basis correspond to physical vector currents, and thus should have vanishing coupling to physical charged 
pions. It is useful to introduce a definition of twist angle incorporating this result: 
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FIG. 1: Illustration of different methods for working at maximal twist. See text for description [method (iv) is discussed in 
sec. [3- The arrows represent the direction one moves to approach the chiral limit (which occurs at /io = 0). The plot represents 
the GSM regime in which m q ~ a. The Aoki regime, in which m q ~ a 2 is represented by the shaded region. It is discussed in 
sec. II VI The angle 8ui is defined to be positive for the situation shown in the Figure. 



where the currents are in the twisted basis, and the superscripts are adjoint SU(2) indices. This angle is denoted 
uj in Ref. but I add the subscript "A" here for clarity. Maximal twist can then be defined by lja — ±i"/2, 
corresponding to 

(OKJO-O, (5) 

Note that the charged pions can be unambiguously created by the local pseudoscalar density (which is invariant under 
the axial transformation that rotates between twisted and physical bases). The condition is equivalent to setting 
the PCAC mass (defined in sec. IIIII below) to zero. 

The two methods apply the condition JHJ) in different ways. In method (i) the condition is enforced for each choice 
of /io by varying mo [13. l23|. This leads to a curve in the mass plane which one follows as /io is reduced. To good 
approximation this curve is found to be a straight line for the masses used in present simulations |23l l24j . This method 
is sometimes called the parity-conservation definition, but I do not use this name as parity is not restored along the 
line it defines at non-zero lattice spacing [and the name could equally well be applied to method (iii) discussed bel ow] . 
Method (ii) defines the critical mass by extrapolating the curve along which the condition (|5jl holds to /io = j24| . 
The resulting simulations are then done holding mo fixed at this value, with /io varying. This method is also that 
proposed in Ref. |r7|. Figure ^ illustrates the difference between the two methods. 6 

The curve defined in method (i) is displayed in Fig. 7 of Ref. |2^| and Fig. 1 of Ref. [24|. (In both figures the 
mo and (jlq axes are interchanged compared to Fig. ^ and in the second figure the vertical axis is 1/ k rather than 
m = 1/(2k) — 4.) These are results from quenched simulations, but I will not rely on the details of these results, 
instead assuming only that the following features are general and continue to hold in unquenched simulations: (1) 
the curves are almost linear for small /i (roughly for no/a < 30 MeV); and (2) they approach the /*o = axis at an 
angle that differs from 7r/2. Indeed this difference (denoted by Scu in Fig.0 is significant: approximately 0.8 radians 
for /? = 5.7 jH, 0.5 for /3 = 5.85 and 0.35 for /3 = 6 H3- 

The first observation I want to make is that the result Slo =/= is a prediction of tm^PT, and that the value of 
5lj is a direct measure of O(a) discretization errors. This follows from results contained in Ref. |13| . In particular, 
method (i) is the same as using the canonical definition of maximal twist suggested in Ref. |13| . It follows from setting 



Both methods (i) and (ii) [as well as method (iii) discussed below] are optimal in the sense discussed in the introduction — they lead to 
automatic O(a) improvement even in the Aoki regime. This property is not important, however, in this section. 
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u>A = 7r/2 in eq. (58) of that work (where, recall, lo a is denoted u>) that the line of maximal twist in method (i) lies 
at an angle 

ujq = 7r/2 - 5uj , 5uj = -(Zm/Z^Sw + 0(a 2 ) , (6) 

relative to the origin where j«o = and m /a = m' c [i.e. ujq is defined as in eq. Here S\y is defined in eq. (J3J 

above, and the factor of Zm/Z^ = Zp/Z$ arises because the predictions of Ref. jljj are given in terms of physical 
masses, rather than the bare masses used in Fig.^ An implication of the result JHJ is that the line defined in method 
(i) extrapolates to a critical mass which includes the 0(a) shift, confirming that it is an optimal choice. It then follows 
that method (ii) corresponds to setting luq — it/2. 

The most noteworthy feature of eq. JHJ) is that the effect is linear in a. The angle 8ui is an example of a quantity 
which is not automatically improved at maximal twist (as discussed more extensively in Ref. |!3jL The result 10 also 
incorporates the result that using different lattice axial currents in the criterion JJJ) leads to different results for Suj. 
For example, using an improved current with ca ^ instead of an unimproved current can be seen to lead to a change 
in Sui proportional to acA- 7 This effect enters through the dependence of Sw on the LEC Wio, for Wio itself depends 
on the choice of current. Note that if one uses non-perturbatively improved Wilson fermions (which sets W = as 
discussed in Ref. and a non-perturbatively improved axial current (which sets Wio — 0), then 8uj is predicted 
to vanish at 0(a). In this case the difference between methods (i) and (ii) appears at next-to- next-to-leading order 
(NNLO) in tm%PT, arising, for example, from a contribution proportional to afi 2 'tp4' i n the Symanzik effective action. 
This latter term gives the critical mass a quadratic dependence on jj,^. The resulting small curvature can in fact be 
seen in the quenched results discussed above, particularly in the results at the largest lattice spacing |24j . 

It should be stressed that a non-zero value for 8u> does not contradict automatic improvement as long as it is of 0(a). 
Contributions to physical quantities linear in a are necessarily also proportional to Suj, so the overall discretization 
error remains of 0(a 2 ). Indeed, as stressed in Refs. there is an intrinsic, irreducible ambiguity in the twist angle: 

different criteria, equally good in the continuum limit, lead to results differing at 0(a) at non- vanishing lattice spacing. 
In fact, one can define maximal twist as lying along any line in the mass plane whose angle satisfies luq = 7r/2 + 0(a), 
where the O(a) offset is arbitrary. Note that for all such lines the quark mass is m q — /x[l + 0(a 2 )], and thus, to NLO 
accuracy, one can use m q — //. 

The above-mentioned curvature in the maximal twist line of method (i) allows one to estimate the size of the 
ambiguity in the resulting critical mass. The NNLO curvature means that, when extrapolating to fj, = from 
fj, ~ O(a), the angle 8ui has an error of 0(a 2 ). This leads to an error in the intercept m' c of size 0(a 3 ), consistent 
with the size predicted by tm^PT and noted in the introduction. 

For practical applications it is important to study Slo as an indicator of the size of discretization errors. One expects 
\6u>\ = a A, with A a scale lying somewhere in the range from Aqcd « 300 MeV to A x = 2V2/ « 1.2 GeV. For the 
(3 = 6 results of Ref. |H (5u> « 0.35 at a" 1 w 2 GeV) one finds A « 0.7 GeV (assuming Z M /Z m « 1). The results 
at the other lattice spacings are roughly consistent with this scale, or, to put it another way, are roughly consistent 
with the expected linear dependence on a. One should not expect more than a rough consistency because there is 
additional implicit logarithmic dependence on a, entering, for example, through the factor Z^/Z mi as well as NNLO 
terms of 0(a 2 ). 

Given the result for 8ui, one can very roughly estimate the size expected for the 0(a 2 ) errors in physical quantities — 
they should be approximately be of relative size (Suj) 2 . This is approximate because different LECs and numerical 
factors enter for each quantity. Nevertheless, if one uses it as a guide, one expects 0(a 2 ) errors of relative size ~ 10% 
for a -1 w 2 GeV. 8 It should be borne in mind, however, that changes to the gluon and fermion actions, as well 
as changing from quenched to unquenched simulations, will change the size of the LECs and thus the discretization 
errors. 

It is interesting to ask what happ ens to the relationship between methods (i) and (ii) for negative fiQ. I n this case, 
it follows from eq. (58) of Ref. 13] with ui = —tt/2 that the figure is symmetric under reflection in the mo axis. In 
fact, this follows from the invariance of the underlying lattice theory under parity transformation combined with a 
change of sign of hq. Thus the line determined by the condition J5j appears to have a cusp at the /io = axis. How 
this extends into the Aoki regime is one of the questions addressed below in sec. II VI 

Knowing the combination of LECs appearing in 8u> allows one to predict other quantities using tmxPT. In particular, 
one of the parity violating axial form factors of the pion [the p\ term in eq. (91) of Ref. [13| at maximal twist and 



7 I thank Andrea Shindler for alerting me to this dependence. 

8 This applies for quantities that do not vanish in the chiral limit, such as f n and mjv- For the pion mass-squared, and in particular for 
the minimum value of the charged pion mass-squared in the case of the first-order transition scenario, or the splitting between charged 
and neutral pion mass-squareds, the 0(a 2 ) effects are additive, and of size a 2 A 4 . 
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for q 2 — > 0] depends only on 8w- Thus, as long as one knows Zp/Z$, one can predict this form factor up to NNLO 
corrections. 

I now turn to my second observation about the relationship between methods for determining maximal twist. A 
third definition of the twist angle is given in Ref. |13|. and denoted u>p. It is based on requiring that the physical scalar 
density does not create the neutral pion from the vacuum. Method (iii) for determining maximal twist is defined by 
setting ojp = ir/2, and is equivalent to requiring that (0|P 3 |7r°) = 0, where P 3 is the neutral pseudoscalar density in 
the twisted basis. In Ref. [13j it was argued that this method would be difficult to implement in practice, since the 
required matrix elements involve disconnected quark contractions. My observation here (which follows from eq. (98) 
of Ref. 01) is that at maximum twist [defined using any of methods (i-iii) ] 

+ 0(a 2 ) = tt/2 + 0(a 2 ) (maximal twist). (7) 

In other words, up to NNLO corrections, methods (ii) and (iii) are the same. Thus it turns out not to be difficult to 
implement method (iii) in the GSM regime. The generalization of this result to the Aoki regime is discussed below. 



III. COMPARING TM X PT EXPRESSIONS TO SIMULATIONS 

In this short section I point out that existing simulation data exhibit qualitative features predicted by tm^PT, and 
that by doing more detailed fits one could over-determine the LECs associated with discretization errors. I illustrate 
these points by showing plots of physical quantities for a choice of LECs roughly consistent with the numerical results 
of Ref. For most of this section I work in the GSM regime, although the plots extend into the Aoki regime, and 
use results from the following section. I assume that an optimal critical mass m! c has been determined, i.e. that one 
knows the critical mass with errors of 0(a 3 ). I then consider the functional dependence of physical quantities on the 
untwisted quark mass ml at fixed values of the twisted quark mass fi. 

As explained in Ref. [2j (and extended into the GSM regime in Refs. E1HD3) 0(a) errors in physical quantities 
can be canceled by "mass averaging", i.e. by suitably averaging the results at (m',/i) and (— to', — /x) (which corre- 
sponds to averaging ±m q in the physical basis used in Ref. 2]). Here the term "suitable" indicates that, for quantities 
with negative R$ parity (defined in Ref. 0), the term at (— m', — fx) must be included with a negative sign 0. Physical 
masses and matrix elements, however, have positive R§ parity. One can change the average to one between (m',/i) 
and (— ml, +/i) by doing, for example, a parity transformation on the second term. This transformation, which I refer 
to as "twisted parity", is equivalent to uj — > — u> in the physical basis. As discussed in Ref. |2j, it does not effect 
physical quantities, but can lead to a sign change in certain correlation functions. The net result is that for physical 
quantities one cancels O(a) errors by averaging over (ml , ^i) and (—to', h) with relative positive signs, while for other 
quantities one must determine the relative sign of the two terms on a case-by-case basis. Two examples used below are 
the PC AC mass and the expectation value (-P 3 ). The former has odd R$ parity, but even twisted parity (at p = 0), 
and so requires a relative minus sign for the two contributions. The latter has odd R§ and twisted parities, and so 
the average has relative positive signs. 

The point I want to make here is that the contributions that are canceled by the mass average described above 
are also of interest, for they provide a measure of discretization errors, and allow tests of tm%PT at NLO. These 
non-continuum parts can be picked out using 

A qln , _ QK,M)-(-) p Q(-m» 

Ab ^>- Q( m >, fX ) + (-) P Q(- m ', fi )> W 

where Q is the quantity of interest, and (— ) p is the product of its R$ and twisted parities. The name refers to 
"antisymmetric", since for most quantities (all those considered here except topcac) AS is the antisymmetric part. 
Using the results of Ref. [13,1 find 

AS(m 2 ± ) = 2AS({0\P*\ir ± ))=AS({Tr\S°\Tr)) = (m'/m q )(26 w -5 w ) + O(a 2 ), 
AS(f A ) = (m'/m q )S w /2 + 0(a 2 ), 
AS(topcac) = (m q /m')S w + 0{a 2 ), (9) 
ASm)) = 0(a 2 ), 
AS(ml -ml±) = 0(a). 

Here /a is the usual pion decay constant, and TOpcac is defined in the following section. I stress that these results 
hold in the GSM regime, but not in the Aoki regime. One sees that tm%PT relates the "antisymmetries" in the 



7 



various quantities. In particular, if Sw has been determined using the angle 8u> as described in the previous section, 
then only one additional parameter (8 W ) is needed to describe all seven quantities in eq. ©■ This assumes that the 
ratio Zp/Zg is known so that bare mass ratios can be converted into those of physical masses. 

One peculiar feature of the results in © is the appearance of m q /m' — 1/coswo m the result for mpcAO By 
contrast, all the other results that are non- vanishing at NLO are proportional to coswoj and thus vanish when ml = 0. 
The divergence at to' = in AS'(topcac) results from the fact that the average is not taken with respect to the 
position at which TOpcac vanishes (which, as seen in the previous section, occurs along a line with loq = tt/2 — Slo 
and not along the line luq = 7r/2). In any case, the result becomes invalid when coswo becomes of O(a), i.e. when one 
enters the Aoki re gim e, so that one does not in fact reach the divergence. 

Figure 5 of Ref. [22 shows plots of ttv^ ± and TOpcac as a function of Too at /Iq = and 0.01. These results are from 
dynamical simulations with a lattice spacing of a -1 ~ 1 GcV. In physical units (and assuming that Zg ~ Zp w 1) 
the twisted quark mass is thus and 10 MeV, while to' lies in the range —100 MeV < to' < 50 MeV. Since 
aAq CD w 90 MeV and <i 2 Aq CD as 30 MeV, the results lie in both Aoki and GSM regimes. Two features of these 
results are relevant here. First, there is a clear lack of symmetry in the results for to 2 ± , and of antisymmetry in 
TOpcac, i-e. the quantities AS(nr^±) and AS (mpcAc) are non- vanishing. Thus it would be interesting to do a more 
detailed fit to determine them. Second, since a significant fraction of the data points lie in the Aoki region, as shown 
by the evidence for a first order phase transition, a comprehensive fit requires a functional form that is valid in both 
GSM and Aoki regimes. The expressions in eq. 10 do not suffice, for they are valid only in the GSM regime. 

The following section is devoted to providing functional forms valid in both regimes. Since the "antisymmetry" 
arises at NLO in the GSM regime, a consistent form must be valid also at NLO in the Aoki regime. In fact, once one 
works at NLO, it may be that the best approach is to simply fit to the full formulae for both positive and negative 
quark masses rather than use the "antisymmetries" directly. 

I close this section by illustrating the forms that result from the joint NLO analysis in GSM and Aoki regimes. I 
use the following parameters, which are roughly chosen to match the behavior seen in Ref. [2^]. The LECs describing 
discretization errors are set to Sw = 5^ — —0.3, 2w' — (0.25 GeV) 2 , and Ws^ — (W3,3 being the additional non- 
trivial LEC needed at NLO in the Aoki region, as discussed in the following section) . In addition I set the continuum 
LECs (the Li of Gasser and Leutwyler) to zero and drop chiral logarithms. Although both these contributions are of 
NLO in the GSM region, and thus should be included for consistency, they do not contribute to the "antisymmetries" . 
Furthermore, they are of NNLO in the Aoki regime. Note that if I define scales by Sw — — aAw, 5 W = —aA w , and 
2w' — a 2 A^,, and use a" 1 = 1 GeV, the parameters I use correspond to Aw = A^ = 300 MeV, and A w > = 500 MeV, 
which are reasonable values for the scales of discretization errors. The choice = is made for simplicity. Taking 
w' > implies that one is in the first-order scenario of Ref. , consistent with the numerical results. The quantity 



\/2w' = 250 MeV is then approximately the minimum pion mass. See Fig. 4(b) for a sketch of the phase diagram. 

Figure [5] shows the forms of to 2 ± for /x fixed at 0, 10, 15 and 25 MeV, and Fig.|3|shows the corresponding plots for 
wpcac- In both figures the left-hand plot shows the results primarily in the GSM regime, while the right-hand plot 
zooms in on the Aoki region. Within the GSM regime, the mass to" (defined in the following section) is equivalent, at 
NLO, to the mass ml used in this section. The prediction of eq. © that the "antisymmetries" should be independent 
of n holds to good approximation for \m"\ > 50 MeV. The values of AS for both quantities asymptote at large |to"| 
to —0.3, roughly the values one obtains from the data of Ref. [22j. 

I discuss the features in the Aoki regime in the next section. I only note here that the first order phase transition 
can be seen in the jumps for fj, = and 10 MeV. The end-point of the transition is for \i = w' /Bo = 12.5 MeV (as 
discussed in the following section) so that the lines with [i = 15 and 25 MeV pass above the transition. 



IV. THE AOKI REGIME AT NLO 



In this section I extend the tm^PT analysis in the Aoki region to NLO. The motivation for making this extension 
is twofold. First, it allows one to see what happens in detail to the phase transition line, and to the lines of maximal 
twist, as one enters the Aoki region. In particular, it clarifies the size and parametric form of the errors that are made 
in different schemes for determining maximal twist. Second, it provides NLO results that can be used consistently in 
the Aoki and GSM regimes. 

The main results of this section are that a NLO analysis in the Aoki regime requires only one non-trivial additional 
LEC, and that this leads to changes in the phase diagram, as summarized in Fig. These changes are small, 
proportional to a 3 as compared to the size of the phase boundary (~ a 2 ), but lead to new effects which might be 
measurable, e.g. curvature in the phase boundary, and a discontinuity in the pion mass across the boundary. I also 
determine the lines of maximal twist, and results for a number of physical pionic quantities, both at NLO. These are 
of theoretical interest, and may be needed to describe the results of numerical simulations. Finally, I describe what 
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FIG. 2: rn 2 7r ± as a function of m" for fixed /i = 0,0.01,0.015,0.025, for the parameter set described in the text. Here m" is 
the NLO quark mass in the Aoki regime, defined in eq. 11211 . For the GSM regime [i.e. most of fig. (a)] m" is equivalent to 
ml at NLO accuracy. All quantities are in units of appropriate powers of GeV. The curves in (a) can be identified from the 
enlargement (b): the lines decrease in thickness as fi increases, with the /i = 0.025 curve dashed. 
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FIG. 3: mpcAC as a function of m" , with parameters and legend as in Fig. [5] 



happens to the phase diagram if the gauge action is tuned so as to reduce the size of the LEC which determines the 
phase structure at LO from ~ a 2 to ~ a 3 . 

In the Aoki regime, LO terms are of size m q ~ a 2 (terms linear in a having been absorbed into the quark mass), 
and NLO contributions are suppressed by a power of a. Thus NLO terms are of size m q a, p 2 a and a 3 . The first two 
are also of NLO in the GSM regime, and thus are already included in the NLO chiral Lagrangian of Ref. [T^- The a 3 
terms are new. Using the properties of SU(2), I find that there are two independent terms: 



/•(» ) 



W3,i r 

p 



Tr^ijTr^E) + h.c. 



3.3 



P 



Tr(itS) 



h.c. 



(10) 



where h.c. indicates hermitian conjugate, and A is the spurion associated with discretization errors, which is set to a 
times the identity matrix. The new LECs, W^ t i and ^¥3,3, are dimensionless. 
Recalling the form of the leading order mass term 



-^Tr( X 't£) + /i. C .. 



(11) 
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where \' — rh' + iT-iji is the mass spurion, we see that the Wa t i term can be absorbed by a redefinition of x' and m'\ 

x'^x" = x' + ^^ = m" + *T 3 £. (12) 

This is equivalent to an 0(a 3 ) shift in m c . Since one does not know m c a priori, this shift does not increase our 
ignorance. All the formulae quoted above, and in Ref. [l^, hold to their stated accuracy with \' replaced by x" (and 
thus M' replaced with M" — y/ rh" 2 + jl 2 and using tanwo = n/m"). 

Thus there is only one non-trivial new term to include in the analysis. Before describing how this changes the LO 
results, I want to comment on redundancy in the LECs. In Ref. [ljj, it was shown that the LEC W\$ could be set to 
zero by a change of variables, as long as other LECs were adjusted accordingly. This implied that, if one kept Wio in 
the analysis, the results would depend only on certain linear combinations of LECs. Extending this analysis to the 
0(a 3 ) terms, I find that the contributions of the two new constants to physical quantities must appear in combinations 
2Wa l3 + Ws,i and W 3) 3 — 32W'W. I have used this to check the results given below, but, to simplify expressions, I 
have set W\q to zero having made these checks. 

Now I turn to the results. My aim is to map out the phase structure at fixed a as a function of m" and \x, and 
to determine NLO expressions for physical quantities throughout this mass plane. Note that, although m" is not 
known a priori, one can use the results given here to determine ml a posteriori. Where appropriate I will include 
contributions which, while of NNLO in the Aoki regime, are of NLO in the GSM regime, and thus are needed for a 
complete, continuous description within the mass plane. 

I begin with the determination of the expectation value, S m , of the pion field £ S SU(2). This is the value which 
minimizes the classical potential, 9 which is here 

-^r = -(c m m" + s m /t)(l + c m 8 w ) - c 2 n w' - c z m w z + 0(a 4 ) , (13) 

r 

where I use 



c m = Tr(S m )/2, and s m = -iTr(£r 3 )/2 , (14) 

and define the useful quantity 

16Q 3 ^3,3 

W 3 = Jg • ( 15 ) 

The range of the variables is — 1 < c m , s m < 1, with the constraint —y/l — < s m < — c^,. For (i ^ 0, the linear 
dependence of V x on s m implies that s m is pushed to one or other end of its allowed range, so that s 2 m + c 2 m = 1. This 
implies that E m lies in the T3 direction. For /t = (the Wilson axis), one can choose the expectation value to lie in 
the T3 direction. Thus from now on I set S m = exp(iLu m T^), so that s m — sinw TO and c m = cosw m . 

In the Aoki region, the LO potential, obtained by setting Sw = W3 — 0, has been analyzed in Refs. p^ . ll5llT^ |. There 
are two scenarios depending on the sign of w' , as shown (including NLO effects) in Fig. For w' < 0, a transition 
line lies along the Wilson axis, within which parity and flavor are spontaneously broken (the Aoki phase [25]), with 
endpoints at rh" — =p2u/. This I call the "Aoki-phase scenario". For w' > 0, there is a first order transition line along 
the n axis with second-order endpoints at jl = ±2w/. This I call the "first-order scenario". 

The extension of this analysis to NLO is straightforward in principle. The extra terms are, by assumption, small, 
and do not change the qualitative features of the two scenarios. They do, however, lead to small changes in the 
positions of the transition lines and the predictions for other quantities. 

I consider first the Aoki-phase scenario (w' < 0). Along the Wilson axis, the —c m m" term in the potential ensures 
that c m is forced to one or other end of its allowed range for large enough \rh"\: c m = sgn(m"). The phase transition 
occurs when the other terms cancel this force, i.e. when dV x /dc m (c m = ±1) = 0. This occurs at 

rh" = =f2u>' — 3ui3 + iw'Sw + 0(a i ) (Aoki-phase endpoints) . (16) 

The first term on the r.h.s. is the LO result, which is of 0{a 2 ). The next two terms are the 0{a?) corrections, which 
lead to a shift in the phase transition line, but do not change its length. Of course, to measure this shift requires 



Loop effects will be shown below to be of NNLO. 
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knowing where the origin is, and I return to this point below. Within the Aoki phase, the identity component c m 
varies continuously from +1 to —1 (in a way which can be determined, if needed, by minimizing V x ). I only note 
here that the expectation value points purely in t 3 direction (c m = 0) when rh" = 0, which is not the mid-point of 
the transition line (unlike at LO). As at LO, the charged pions are massless in the Aoki phase, with the neutral pion 
also massless at the end points. 

The changes at NLO are more significant in the first-order scenario. There is no symmetry which requires the 
transition to lie parallel to the /i axis once NLO terms are included. To begin I consider the position of the transition 
line on the Wilson axis. The transition occurs when the two minima at c m = ±1 have equal energy. Because the w 3 
term in V x is an odd function of c m , this occurs not at rh" = (the LO result) but at rh" — — w 3 . For rh" greater 
than (less than) this value, c m = +1 (—1). 

This analysis is easily generalized off the Wilson axis. At LO, the first-order transition is present at rh" = if 
\£l\ < 2w' . The LO potential has equal minima at c m = ±c^ = ±yl — (s^,) 2 , both having s m = = fi/(2w'). 
The two minima coalesce into second-order cndpoints at a = ±2w'. At NLO, the terms odd in c m give opposite 
contributions at the two minima, and one must adjust rh" to make their energies equal. To the order I work, rh" can 
be determined by enforcing V x (c m ~ = V x (c m = — c° n ). From this it follows that the transition is shifted to 

rh" = -w 3 + (w 3 - 2w'S w ) ( Jh) + 0(a 5 ) (1st order line: < 2w') . (17) 

\2w' J 

The transition line is thus a quadratic function of /i, which maintains the symmetry under fi — > — /i while being 
smooth at /1 = 0. The value of fi at the end-points is unchanged from LO, aside from NNLO corrections of 0(a 4 ) 
discussed below. For a given jl, however, the error in rh" turns out to be of (3(a 5 ), rather than the naively expected 
0(a i ). Finally I note that the NLO terms shift the angle uj m at the transition line from the LO value by 

Sco m = w m - < = (S w - w 3 /w')(u-/[2w']) + 0{a 2 ) . (18) 

This 0(a) shift is in the same direction on both sides of the transition, so the resulting values of c m on the two sides 
differ in magnitude: c m = ±cj^ — (Suj m ) s^. This leads to the discontinuity in m 2 at the transition, as shown in 
Fig. |2(b)| and discussed below. 

To understand the size of the uncertainty in the position of the end-points it is useful to display the form of the 
potential there. The considerations are identical for both end-points, and I choose to consider that at positive /z. 
Setting fi = 2w' and rh" = —2w'8w [from eq. (|17|) ]. one has 

V w'c 4 

= + {w'Sw - w 3 )c 3 m + 2w'8'^ v c 2 m (1st order endpoint). (19) 

I have dropped terms of higher order than quartic in c m , as they would be subleading in the following discussion. At 
LO, only the first term on the r.h.s. contributes. It has a minimum at c m = 0, and corresponds to the end-point since 
the curvature vanishes. At NLO, the value of rh" has been chosen by the criterion given above so that c m = remains 
a stationary point of the potential. c m = does not, however, correspond to an end-point, for there is a non- vanishing 
quadratic term at c m = (with coefficient cx a 4 ). Furthermore, for some values of the LECs, the cubic term leads to 
a second stationary point with c m ~ a. To move to the actual end-point, one must shift ji and, in general, rh" . For 
example, the quadratic term can be canceled by a shift of size Sfi ~ a 4 , because 5V X cx Sfi(l — c^/2 + ...). It can be 
shown that a similarly sized shift in ft, as well as a shift Srh" ~ a 5 , will also remove the second minimum, if present. 

A NLO calculation does not, however, control the size of these shifts. This is because there are neglected higher 
order contributions to the potential which are of the same size as those in eq. (I19H when c m ~ a. For example, there 
will be NNLO and higher order contributions to the classical potential of the form SV X ~ a^<? m + a 5 c m , depending on 
additional unknown LECs. These, like the terms in l|19ll . are of 0(a 6 ) when c m ~ 0(a), and thus compete with NLO 
terms. In particular, the term linear in c m leads to an unknown shift Srh" ~ a 5 , and the quadratic term to a shift 
6£l ~ a 4 . 

In addition, loop contributions to the effective potential enter at the same order. Flavor-conserving four-pion 
vertices of LO (and thus proportional to p 2 or m q ) give rise to quadratic terms in the potential of the form 5V X ~ 
c 2 m rn\ ln(m w ) ~ a 4 c^. The latter equality follows because, for charged pions, m 2 ~ a 2 at the end-points. There 
are also flavor- violating four-pion vertices of the form a 2 s 2 n W'ir ■ ifir 2 (see Ref. 13]), and these also lead to SV X ~ 
c 2 n a 2 m 2 lnm T ~ o^c 2 ^. Finally, there are tadpole diagrams involving three-pion vertices. These vertices have been 
given in Ref. in the GSM regime. In the Aoki regime they become 

/4* = -6 W 7T 3 d,7f ■ d.TT + Sm7% " ~ Cm *V 3 7T • 7? . (20) 
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FIG. 4: Illustration of phase structure and lines of maximal twist at NLO in the Aoki regime for both Aoki-phase and first-order 
scenarios. The solid lines indicates first-order transition lines, with the solid circles being second-order end-points. Both masses 
m" and /i range up to values of 0(a 2 ) in this figure, with other features (the asymmetry in the phase line in the Aoki-phase, 
the offset of the ujp = ±7r/2 lines, and the horizontal extent of the phase-line in the first-order scenario) being of 0(a 3 ). See 
the text for more details. The angle 8u> is the same as that in Fig. 



Noting that the coefficient of the second term is of 0(a 3 ) at the end-points, I find that the contributions to the 
potential are of the form 

5V* ad ~ c m {am\ lnm^ + a 3 ml lnm T ) — a 5 c m , (21) 

and thus of the same size as the higher order terms in the classical potential for c m ~ a. Clearly extending the 
calculation beyond NLO will be challenging. 10 

A corollary of the previous discussion is that, close to the end point, the NLO calculation does not control the shift 
Su> m , since higher order terms contribute with the same size as LO and NLO terms. For the result (|18fl to be valid, the 
distances from an end-point must satisfy 5(1 3> a 4 or 5m" 3> a 5 . Furthermore, as one approaches the end-points the 
rate of convergence of the chiral expansion slows. For example, if one moves a distance 5m" ~ a 3 from the end-point, 
then one can show that c m ~ a 1 / 3 . It follows that each higher order is suppressed by a relative factor of ~ a 2 / 3 , 
rather than ~ a. 

The results for the phase structure are summarized in Fig. 0] The circles around the end-points emphasize that 
the chiral expansion breaks down in their vicinity, as just discussed. Note that, since the origin of m" is not known, 
in the Aoki-phase scenario one can, at this stage, only determine w' from the length of the phase region, while in the 
first-order scenario one can determine both w 1 (from the length in the fi direction), and W3 — 2w'5w (from the extent 
of the phase line in the m" direction). 

I now turn to results for other physical quantities. To do so, one needs to know the value of io m for a general 



The breakdown of power-counting near the end-points has been noted in the Aoki-phase scenario in Ref. 12(1 . 
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position in the mass plane. This is determined by solving 

= 7^-7-^ = s m m" - c m fi + 2c m s m w' + (2c m s m m" + [s 2 n - c 2 n ]fi)5w + 3c 2 n s m w 3 + 0(a 4 ) , (22) 
j dui m 

and choosing the solution with minimal energy. In practice, for general values of the parameters, and away from the 
phase boundaries, one can solve the LO equation (the first three terms on the r.h.s.) and use this solution as the 
starting estimate in a numerical solution to the full NLO equation. 

I first consider the lines along which uia = ±7r/2, i.e. the lines of maximal twist as defined by method (i) of seclTfl 
The advantage of choosing variables such that Wio — is that uja — ^>m + 0(a 2 ) [l^- Thus the task is to solve for 
ui m — ±7r/2, so that c m — and s m = ±1- For these values, eq. (1221) simplifies considerably, leading to 

m" = [t5 w + 0{a 2 )] ft {lo a = ±tt/2). (23) 

In the Aoki-phase scenario, this result holds for all fi, while in the first-order scenario it holds for \fi\ down to a 
minimum value of 2w' , at which point the uja = ±7r/2 lines run into the second-order end-points. 11 For smaller fi one 
can show that the solution l|23|) does not have minimum energy. The situation is shown in Fig. 0] 

The result eq. (12; jfl answers the question (posed in sec. [HJ| of what happens to the lines of maximal twist when 
one enters the Aoki region. The answer is that they remain straight [the angle given in eq. 10 agrees with that in 
eq. (I23[) to the stated accuracies] , and extrapolate to the origin in the mass plane as long as one uses the variable to" 
rather than to'. In other words, the 0(a 3 ) uncertainty in ml in the GSM regime (discussed in sec.|nj| is resolved by 
extending the calculation into the Aoki regime. 

Note that in both phase-transition scenarios one can determine 5\y, w' and W3 from knowledge of the position of 
the phase transition lines and of the line on which uja = tt/2. In both cases, Sw comes from the slope (assuming 
Zs/Zp is known), and w' from the extent of the phase transition line. In the Aoki-phase scenario, the phase line is 
asymmetric with respect to the origin, and from this one can determine W3 using eq. (|16|) . In the first-order scenario, 
one can use the curvature of the phase boundary or, equivalently, its intercept with the Wilson axis, to determine W3, 
using eq. l(T7l) . 

I next determine the lines along which uip — ±tt/2. It is shown in Ref. [l3l ] that ujp = ui m + s m 5w + 0(a 2 ) in both 
GSM and Aoki regimes. Thus one needs to solve eq. I|22|l for the lines along which ui m = ±tt/2^fSw- At the accuracy 
I work this means that s m — ±1 and c m = S\y. The result is very simple: 

m" = -2w'S w + 0(a 4 ) {u> P = ±tt/2). (24) 

Recall that in the GSM regime, the uip = ±ir/2 line lies on the /i axis for any suitably accurate definition of the critical 
mass. In particular, one could extrapolate the uja — ±7r/2 lines to fj, = and use this to define m c [i.e. method (ii) 
of sec. m. The result (|24H shows that, with the greater resolution provided by working at NLO in the Aoki regime, 
there is an 0(a ) offset between the definitions of m c determined from uja and uip. In other words, within the Aoki 
regime, methods (ii) and (iii) differ at NLO. 

The situation is illustrated in Fig. ^ Note that, in both scenarios, the position of the uip — tt/2 line is a prediction 
of tm^PT, since it depends on the combination Sww 1 that can be determined as discussed above. In the Aoki-phase 
scenario, the lines of uia = tt/2 and uip = tt/2 do not cross for either sign of 5w- In the first-order scenario, the lines 
do meet, and they do so at the end-points of the phase transition line. In fact, since the ui m = tt/2 solution ends at 
these points, so do the uip = tt/2 lines. 

I have focused on the lines of maximal twist since these are of greatest practical interest. I note, however, that one 
can predict the values of uia and uip to NLO accuracy using the formulae given above throughout the Aoki and GSM 
regimes. In particular, the jump in |c m | as one crosses the first-order line [see eq. I|18fl and subsequent discussion] 
applies also to \ca\, since uia — ui rn . 

I now turn to pion properties. I find the charged pion masses to be 

ml± = (c m m" + s m /t)(l + c m [2S w - + 2c 2 m w' + cf„(3w 3 - 25^w') + M" 2 terms + 0(a 4 ) . (25) 

The "M" 2 terms" are as in the continuum, with 2B^m q replaced by M" = \J ml' 2 + jl 2 . They include chiral logarithms, 
and are given, for example, in Ref. [2]]]. Within the Aoki regime they are of 0(a 4 ) and thus of NNLO, but they are of 



11 There is an apparent inconsistency here: the end-points both have uim = uia = ±tt/2 and, according to eq. 1181 . u> m = t /2 + Sui m - This 
is resolved by the fact that the O(a) corrections to uim are not controlled in the present calculation near the end-points. 
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NLO in the GSM regime and must be included there. The new NLO terms in the result (|25|l are those proportional to 
a 3 c^. Given a set of LECs, one can obtain the pion mass at the desired position in the mass plane by first determining 
co m using eq. I122II and then substituting in eq. H25J1 . This is how the plots in sees. IIIII and Ivl are made. 

Although the general formula is complicated, it simplifies at maximal twist. In particular, if one works along the 
lua = Lo m = tt/2 line, i.e. using method (i), then all the NLO corrections vanish: 

ml ± =(i + M" 2 terms + 0(a 4 ) {lo a = tt/2). (26) 

Thus not only are 0(a/i) terms absent, but also 0(a 2 ) and 0(a 3 ) terms vanish. The result thus takes the continuum 
form (with m q = fi, and up to 0(a 4 ) corrections) all the way until the end of the loa = tt/2 line (i.e. until /i = for 
the Aoki- phase scenario, and until /j, = 2w' for the first-order scenario). The same result holds to the stated accuracy 
with both methods (ii) and (iii) applied in the Aoki regime, aside from one caveat. On the ujp — tt/2 line, c m ~ a and 
to" ~ a 3 , so the terms in eq. Ij25(l which differ from the continuum result remain of 0(a 4 ) in the Aoki regime. The 
same is true on the m" = line in the Aoki-phasc scenario. The caveat is that, on the to" = line in the first-order 
scenario, u m starts to differ by 0(1) from tt/2 when /i approaches the end-point. Thus one can only use method (ii) 
in this case down to fi > 2w' . 

Along the Wilson axis, by contrast, there are corrections of 0(afi), 0(a 2 ) and 0(a 3 ). These lead to a complicated 
behavior, including a discontinuity at the first-order boundary: 

m 2 ± = 2w ± 2(w 3 — S^rW 1 ) + 0(a i ) (first-order discontinuity, /Lt = 0). (27) 



The behavior can be seen in Fig. 2(b) which also illustrates the curvature of the transition line because the boundary 
is at non-zero to" for /j, ^ 0. 

One check on the result eq. (|25|) is that, in the Aoki-phase scenario, the charged pion mass vanishes at the end-points, 
and remains zero (within the 0(a 4 ) accuracy) inside the Aoki phase. 
The NLO result for the flavor-breaking splitting Am 2 — m^ — m 2 ± is 

Ami = -2s 2 m w' + 2s 2 m c m (w'[5ft + 25 w ] - 3w 3 ) + 0(a 4 ) . (28) 

Note that the a 3 term vanishes both on the Wilson axis and for io m = ±tt/2. In fact, to NLO accuracy, Am 2 — —2w' 
along both the loa = ±7r/2 and lop = ±tt/2 lines, i.e. for both methods (i) and (iii), as well as for method (ii) with 
the same caveat as discussed for the charged pion mass above. 

The a 3 terms do not contribute to the pion decay constants f a and fp, nor to the scalar form factor of the pion or 
the condensates. Thus the expressions for these quantities given in Ref. remain valid and I do not repeat them 
here. I do stress, however, that at maximal twist, using any of the definitions, all these quantities agree with their 
continuum forms up to NNLO. 

The PCAC mass does, however, receive contributions from 0(a 3 ) terms. The definition that is useful for tm^PT 

is 

_ (Ol^k 1 ) 

WPCAC = o/nlnll U ' ( 29 ) 



This leads to the general expression [ij 



tan loa tan lo, 
which is useful except on the Wilson axis, or to the result 



m PC Ac = — ^ = -^[l + 0(a 2 )], (30) 



-pcac = 5=^1 [1 + 0(a 2 )} = C -^- [1 + c m (S w - M + 0{M") + 0(a 2 )] . (31) 
ZJp its 

These formulae are used to make the plots in Fig. |21 Note that the jump in topcac a t the phase-boundary is 
asymmetric, i.e. the minimum value is different on the two sides of the transition. 



I close this section with a brief discussion of what happens if, by adjusting the gluon and quark actions, one were 
able to reduce the size of w' from 0(a 2 ) to 0(a 3 ). This would be advantageous from a practical viewpoint since the 
size of the region in which lattice artifacts have an important influence on vacuum alignment (the Aoki regime) would 
be reduced. Furthermore, since w' is the only 0(a 2 ) term, discretization errors in general would be reduced. 

The discussion of phase structure given earlier was predicated on the 0(a 3 ) terms being small perturbations to 
the w' term. If this is not the case, then qualitative changes are possible. A general analysis for w' ~ 0(a 3 ) is 
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straightforward in principle (although one must work only at LO, since 0(a A ) terms are not controlled) but I have 
only worked out in detail the simplest case of w' — 0. The vacuum is then determined by competition between cubic 
and linear functions of c m . I find only one scenario: as one moves along the /i = axis there is a second-order 
cndpoint, followed by a region of Aoki-phase with spontaneous flavor-parity breaking of length ~ a 3 , and then by a 
first order transition. This latter transition presumably extends in the ±/i directions for a distance of ~ a 3 , ending 
at second-order endpoints. In other words, the usual two scenarios are merged into one having features of both. The 
remnant of the presence of two scenarios for w' ~ 0(a 2 ) is that the relative positions of the first-order transition and 
Aoki-phase segment depends on the sign of W3. Presumably, as w' increases in magnitude, one or other feature of 
this picture will reduce and then disappear. For example, if w' is positive the Aoki-phase segment will reduce in size 
and then disappear, while the length of the first-order transition will increase. 



V. THE ABSENCE OF INFRA-RED DIVERGENCES 

In this section I comment on the work of Ref . [l^ ■ Based on a general analysis using the Symanzik expansion , 
Ref. |17| finds apparently infrared (IR) divergent discretization errors in physical quantities at maximal twist which 
are at worst of the form (a/m 2 ) 2k , with k > 1 an integer. These divergences are interpreted as indicating the onset of 
large discretization errors as m 2 . is reduced, and the breakdown of the Symanzik expansion when m 2 ~ m q < a. This 
result is for a non-optimal choice of m c having an error of 0(a). Using an optimal choice (such as those discussed in 
sec.|nj Ref. 01 finds the divergences to be weakened, such that the leading IR divergent contribution is proportional 
to a 4 /m 2 .. This leads to the conclusion that it is possible, with an optimal choice of m c , to work with quark masses 
as small as m q > a 2 , although for smaller masses, i.e. in the Aoki regime, discretization errors become large. 

The conclusion that one must use an optimal choice of m c to allow one to work in the GSM regime is in agreement 
with that obtained previously using tm X PT inRefs. [HE^. Thus there is no dispute over how to proceed practically, 
and indeed all recent simulations use an optimal choice for m c [2^, 0j H3- My observations here concern the 
interpretation of the apparent IR divergences, their implications for the validity of the Symanzik expansion at small 
quark masses, and the theoretical status of calculations using tm^PT. 



A. Non-optimal critical mass 

I consider first a critical mass having an error of O(a), since in this case tnr\;PT at LO provides simple examples of 
the fate of the IR divergences. What one learns in this case can then be generalized to the situation with an optimal 
m c . One superficial complication when comparing the results of Ref. [lj| to those of tm%PT is that the former are 
given in the physical basis while the latter are in the twisted basis. Thus I begin by briefly summarizing the approach 
and results of Ref. [17| using the language of the twisted basis. Their first step is to determine the form of the 
Symanzik local effective Lagrangian corresponding to the lattice theory under study. The target continuum theory is 
given by the operators of dimension 4, 

£4=^^ + ^757-3^- (32) 

where my use of the twisted basis shows up in the form of the mass term. Note that fj, is the physical quark mass 
(called m q in Ref. [l7|h The dominant discretization errors come from the dimension 5 contribution, which, in the 
twisted basis, is 

a£ 5 ~ a-ipio^F^ijj + aAQ CD ^ + ay 2 ^ ■ (33) 

Here the notation is schematic, showing only the form and order of magnitude of each term, but omitting the unknown 
coefficient of O(l) which multiplies each term, and which can have a logarithmic dependence on a. The final term in 
£5, being proportional to /1 2 , does not contribute to the leading apparent IR divergences, and can be dropped from 
the subsequent discussion. 

At this stage the analysis can be seen to be identical to that used in the tm^PT approach, as described, for example, 
in Ref. Il6l. There is, however, one subtlety in this comparison. The result in eq. (|33fl apparently differs from that 
of Ref. |16| by the absence of the aAQ CD ?/5-0 term in the latter. This can be traced to the definition of critical mass 
used in Ref. |l6j |: the quantity denoted m c is defined to include all perturbative and non-perturbative contributions to 
the untwisted mass term. In other words, in the approach of Ref. |l6(, the aAQ CD ?/j^> term in £5 is moved into £4 and 
absorbed by an 0(a) shift in the definition of the untwisted quark mass. Conversely, maximal twist in the definition 
used by Ref. |17| corresponds to having an untwisted mass m ~ aAq CD in the £4 of Ref. |l6j |. and this is naturally 
moved to £5. In fact, in Ref. |l6j | there is an additional term in £5 of the form am 2 ?/'?/': but this is now seen to be of 
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size a 3 Aq CD ?/5^, and thus resides in €7. An additional operator in Ref. 0] of the form mipftip — > Aqcd^PP ip can be 
absorbed by a redefinition of the fermion field. 

The tmxPT analysis and that of Ref. [l3| now diverge. In the former it is noted that the operators in £5 transform 
under SU(2) chiral symmetries like mass terms, and thus it is known from the methods of chiral perturbation theory 
how to systematically incorporate their effects into the chiral Lagrangian describing the low-energy physics of QCD. 
At LO, the transcription is 

a£ 5 — ► -a^Tr(£ + E f ) , (34) 

where, as above, a = 2Woa, with Wq ~ Aq CD an unknown LEC. The r.h.s. of eq. (|34[1 has exactly the same form as 
an untwisted mass term, with a mass of size aAq CD , in the chiral Lagrangian. 12 Thus if one studies the theory as 
a function of fj,, one is working along a line of the type denoted "method (iv)" in Fig. ^ In tmxPT, the effects of 
£5 can be completely determined at LO by simply incorporating into the target continuum theory of £4 a n untwisted 
mass 5m = a/(2Bo). The resulting theory has no IR divergences (since the total quark mass m q = a//i 2 + {5m) 2 
does not vanish) and one can obtain expressions for physical quantities valid for all /i. In particular, the remaining 
contributions from £5 [i.e. those of NLO in tm^PT for which the transcription into the effective chiral theory is more 
complicated than that in eq. 1|34[) ]. as well as those from £$ and higher order terms in the Symanzik expansion, lead to 
corrections suppressed by the expected powers of csAqcd El • This holds true also for the properties of particles 
other than pions |28j . 

I return now to the analysis of Ref. [l7j . This is based on the fact that £5 (and £7, etc.) creates a neutral pion 
from the vacuum (recall that tptp m the twisted basis corresponds to iV' p hys75' T "3'!/'phys m the physical basis), as well 
as giving rise to vertices among odd numbers of such pions, and to matrix elements in which other hadrons couple to 
odd numbers of neutral pions. The resulting neutral pion propagators at zero four- momentum give factors of \jm\. 
A detailed accounting of these contributions leads to the conclusions given above concerning the presence and form of 
IR divergences. In particular, if m 2 ~ a, then all the leading IR divergences proportional to (a/m 2 ) 2fc are of the same 
size, suggesting a breakdown in the Symanzik expansion. Note, however, that these leading IR divergences result only 
from multiple insertions of £5, and do not involve higher order contributions to the Symanzik effective Lagrangian 
(£q, £7 etc.). This is an indication that a summation of all apparent IR divergences may be possible. 

Indeed, my observation here is that tm^PT provides such a summation. This is accomplished, as outlined above, 
simply by absorbing the additional 0(a) term into the original mass term in the chiral Lagrangian. Note that since 
5m ~ a one is necessarily in the GSM regime. The shift in the mass causes a change in the vacuum: at LO, the 
twist angle of the condensate becomes uj m = tan _1 (/i/i5m) rather than ir/2. The apparent IR divergences are then 
seen as an indication that one is expanding around the wrong vacuum. Of course, it would be possible to force the 
vacuum in tm^PT to remain at lo™ = 7r/2, and treat the term in eq. I|34[l as a perturbation. This would reproduce 
the diagrammatic analysis of Ref. |l7j . In any case, the main point I want to make is that the summation provided 
by tm%PT shows that there is no breakdown of the Symanzik expansion, in the sense that it is legitimate to treat 
subsequent terms as having progressively smaller contributions. The subtlety here is that one must treat £5 non- 
perturbatively, so that the manifest powers of a associated with £5 are lost. But once this is done, as is possible with 
tmxPT, the manifest powers of a from subsequent terms (£q etc.) are retained. 

Let me illustrate these general observations with two concrete examples. It is shown in Ref. |17| that the most 
divergent quantities are the overall correlation functions, in contrast to particle energies, where the divergence is 
ameliorated by having one less factor of \jm\. Thus I consider the two-point function of the "physical" axial current, 
^i,phy S ; a t large Euclidean separation, so that the charged pion contribution dominates. In the approach of Ref. [l7j . 
one proceeds as though one is at maximal twist, and thus the physical axial current is, in the twisted basis, given 
purely by the vector current: A 1 ^ phys = V*. Thus the correlation function of interest is 

(4, phys (z)< phys (0)) = (V*(x)V*(0)) . (35) 

Here, for simplicity, I am assuming that both lattice axial and vector currents are multiplied by appropriate Z-factors 
so that they are correctly normalized. Using the results of Ref. ^3], the VV correlator can be evaluated in tm^PT. 
It is sufficient for illustrative purposes to work to LO. Then one finds 

(V*(x)V*(0)) = sin 2 u; m (Al(x)Al(0)) + non-pion contributions , (36) 



12 At this point the choice of whether the aAQ CD ^^> term was placed in £4 or £5 becomes irrelevant, since the total coefficient of Tr(E + E^) 
is the same in both cases. 
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fA/f, m" = . 005, .015, .025 pionmass squared, ra" = . 005, . 015, . 025 




(a) f A /f (b) ml ± 



FIG. 5: Examples of quark-mass dependence when using a critical mass with an error of 0(a) [method (iv) of Fig. 0: (a) 
f A /f and (b) ml ± (in GeV 2 ) versus (in GeV) The parameter-set is described in the text, and corresponds to the first-order 
transition scenario, m" is fixed to 8m, with values 0.005, 0.015 and 0.025 GeV. These are distinguished by the line width, 
which increases with the magnitude of Sm. The dashed line shows the continuum result. 



where A is the actual LO physical axial current appropriate to a vacuum oriented such that tancj m = /j,/Sm. The 
overall factor in front of the correlator enters directly in the decay constant (squared), and so one finds, at LO: 

f A (method (iv)) . [i I ( 5m\ 2 

= smu) m = = 1 - - — +.... (37) 

/ vV + {5m) 2 2 \ /i / 

Since, at LO, m 2 oc fx, the corrections to unity have precisely the IR divergent forms deduced by Ref. ^t|- As can be 
seen, however, they sum up into a simple, predictable form. 13 

The second example is the charged pion mass. At LO in tm^PT the result is simple: 

1 (5m) 



m 2 ± (method (iv)) = 2B Q ^fi 2 + (Sm) 2 = £ \l + + • ■ ■ J • (38) 

Thus the corrections to the squared charged pion mass begin at order (dm) 2 / \i, one power of 1/n less infrared divergent 
than the corrections to Ja- This (as well as the form of the higher order terms) is as predicted by Ref. ^tJ- Once 
again, however, the apparently IR divergent terms are summed up straightforwardly. 

Before turning to the issue of apparent IR divergences when using an optimal m' c , I think it useful to show examples 
of the expected behavior of the quantities just discussed for various choices of 5m. To do so I use the full NLO results, 
valid in both GSM and Aoki regimes, rather than the LO expressions for the GSM regime given in eqs. I|37ll and 
This gives a more realistic view of the expected forms. Figures El and El show the results of tm^PT for /a// and rn?±. 
The parameters are as for Figs. Eland El (5w — <fe = —0.3, 2w' = (0.25 GeV) 2 , W3 = 0, Li = and chiral logarithms 
dropped), except that I have added a non- vanishing continuum analytic part to /a (setting L4 + L5/2 = 0.007) to 
make the slope versus quark mass more realistic, and that the sign of w' is changed for Fig. E| so as to show the 
behavior for the Aoki-phase scenario. For my choice of parameters the end-points in the first-order scenario (Fig. EJ) 
lie at (m", fj) « (4, ±12.5) MeV while the Aoki-phase (Fig. El runs from m" = -5 MeV to +20 MeV. This substantial 
asymmetry is due to the relatively large size of 5w- 

In all the plots the dashed lines are the continuum results for the parameters I have chosen, assuming that one is 
at maximal twist. The solid lines illustrate the dependence on the twisted mass /i for three choices of the untwisted 



mass 5m. I use 5m = 25, 15 and 5 MeV, except in Fig. 6(b) where I use —25, —15 and —5 MeV. These values are 
chosen to bracket the GSM and Aoki regimes (recall that a z Aq CD = 7.5 MeV for a^ 1 = 2 GeV). The largest value of 
\5m\ (giving the curves with the thickest lines) is most representative of method (iv), i.e. it has 5m ~ a for typical 
choices of a. The curves for \5m\ = 25 MeV bend away from the continuum form for [i < 5m, as expected from 



At NLO in the GSM regime, the result from method (iv) can be shown to give sinu m times the NLO result obtained in Ref. 1 1. 'I . It is 
this NLO form which is plotted below in Fig. |5] 
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fA/f, m" = . 005, . 015, . 025 pionmass squared, -m" = . 05 , . 1 5 , . 02 5 




(a) f A /f (b) ml ± 



FIG. 6: As in Fig. ^except that the sign of w' is changed so that one is in the Aoki-phase scenario. In addition, the values of 
8m axe chosen to be negative in (b). 

eqs. (|37H and (|38|l . The NLO corrections are significant, however, as can be seen from the fact that the results are 
different for first-order and Aoki-phase scenarios. In other words, the results depend on the choice of w' , which is a 
NLO parameter in the GSM regime, rather than just on \5m\ and /it, as predicted by the LO results. 

I include the results for \5m\ = 5 MeV to illustrate what happens if one enters into the Aoki regime. These do not 
correspond to method (iv) as defined above, and the expressions l|57|) and are not valid when /j, « \Sm\. This can 
be seen clearly from the 5m — 5 MeV curves for the first-order scenario in Fig. [5] they are not smooth, and there 
is a very large difference between them and those for the Aoki-phase scenario. I choose negative values for 5m in 
Fig. so that only the 5m = — 5 MeV result corresponds to running into the Aoki phase — indeed this curve shows 
what happens when one runs into an end-point of the Aoki phase, while those for 5m = — 15 and —25 MeV avoid the 
phase boundary. The deviations from the continuum result are also larger for negative than for positive 5m (another 
indication that the LO expressions are insufficient). 

The intermediate choice |<5m| ~ 15 MeV lies at, or close to, the boundary of the Aoki regime, depending on the 
value of a, and I include it for completeness. 

Although the discussion in this section is mainly theoretical, it is interesting to know whether method (iv) has been 
used in practice, i.e. whether the error in m c is of O(a) or smaller. This is not entirely clear. The traditional method 
for determining m c has been to extrapolate m\ to zero along the Wilson axis with a simple fit function starting from 
relatively large quark masses. There are two issues that arise in evaluating this "m^ = 0" method. First, even if 
one did a perfect extrapolation, the result would differ by 0(a 2 ) from an optimal critical mass Second, the 

long extrapolation introduces a systematic error which, while not parametrically of 0(a), may be numerically of 0(a) 
in a given simulation. For example, in the case where there is an Aoki phase (as appears to be true for quenched 
simulations with the Wilson gauge action), the extrapolation might miss the actual end-point by an amount of O(a). 
If so, one would be using what I am cal ling method (iv) . 

The quenched studies of Refs. |23l Y2A 1271 129| compare the results at maximal twist using the "to^ = 0" definition of 
the critical mass to those obtained with an optimal choice (using methods (i) or (ii), or, in the case of Ref. [2j|, results 
from overlap fermions). The results using the former definition are found to bend away from those with an optimal 
m c for [i < a. In light of the previous discussion, however, it is unclear whether the "m^ = 0" definition corresponds 
to method (iv), and thus one does not know which of the curves in Figs. |3] or best illustrate the expected behavior. 
For example, in the (3 — 6 data of Ref. (2^, the value of 5m is found to be approximately 0.01(Z m /a) ps 20 MeV. 
Given that aAg CD ps 45 MeV and a 2 AQ CD w 7 MeV for this simulation, it is unclear whether 5m should be treated as 
of 0(a) or 0(a 2 ). From a practical viewpoint, however, this is not important. One should attempt to fit the results 
to the NLO tm%PT forms keeping 5m as a free parameter. 

B. Optimal critical mass 

I now return to the case of an optimal m c . The discussion follows a similar line to that for a non-optimal m c , but 
in this case there are no simple analytic forms to illustrate the summation. 

In the analysis of Ref. [l7|], the apparent IR divergences are suppressed because, by construction, the single pion 
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matrix elements of a£$ are reduced from ~ a to ~ afi, a 3 (and are thus the same size as the matrix element of a 3 £? if 
fj, ~ a 2 ). The dominant terms are of the form a 2 (a 2 /m 2 ) k , with k > 1, and so are all equally important when m 2 ~ a 2 , 
at which point they are all ~ a 2 . These terms apparently preclude one from working in the Aoki regime, for although 
the errors are parametrically small (proportional to a 2 ) they have a complicated and non-uniform dependence on the 
lattice spacing and quark mass. Furthermore, the Symanzik expansion is apparently breaking down. 

The tm%PT analysis shows, however, that, as above, the IR divergences indicate expansion about the wrong vacuum 
and can be summed by implementing a non-perturbative shift to the correct vacuum. As noted in Ref. [T^ , the leading 
IR divergences result from multiple insertions of £§, with only one creation of a neutral pion by the suppressed vertex 
£5,7. Correspondingly, in the tm^PT analysis the change of vacuum is that caused by the a 2 coefficient w'. This is 
exactly the analysis of the phase structure in the Aoki regime carried out previously in Refs. m> 0, 0, , and 
extended here to NLO. The results for Ja/J and m 2 are exemplified by the \5m\ = 5 MeV curves in Figs. EJandEJ 

The conclusions I draw are as follows. First, there is no obstacle, in general, to simulating in the Aoki regime. 
The tm%PT analysis provides functional forms that allow one to predict and fit the behavior of physical quantities in 
this regime. There is, in other words, no breakdown of theoretical understanding in this region, and in particular no 
breakdown of the Symanzik expansion, in the sense that once one includes the leading effects of £q non-perturbatively, 
the subleading contributions from these operators, and the contributions of higher order operators (including now ^5), 
are systematically suppressed by the expected powers of oAqcd- In fact, if one works at maximal twist, defined using 
any of methods (i-iii) (with the caveat that one must stay above the end-point in the first-order scenario), then the 
contributions of £q remain perturbative also, along with their manifest powers of a. This is because, when using any 
of these methods, the condensate points in direction of maximal twist up to small corrections, uj m = tt/2 + 0(a). 

Second, although the apparent IR divergences are summed up by tm^PT, their "residue" is a complicated depen- 
dence of physical quantities on a and fi. Third, and most important from a practical point of view, if one works 
at maximal twist using any of methods (i-iii), then the continuum extrapolations are not complicated. Indeed, as 
already noted in the previous section, there are no a 2 (or a 3 ) corrections to the charged pion masses in this case. 
Thus in the Aoki-phase scenario, there is no obstacle to working all the way down to /1 = 0, although in the first-order 
scenario one must stop at \fi\ > 2w'. Nevertheless, the overall point is that one does not need to impose the constraint 

m g ~ n> a 2 mm 

There is an exception to the claim that tm^PT controls the IR divergences. This is in the vicinity of the end-points 
of the phase-boundary. Here, as discussed in the previous section, the power counting of tnr\RT breaks down, and 
higher order terms, including loops, are needed to determine the vacuum. Thus, in the first order scenario, one should 
work at maximal twist only until one is a distance Sfi ~ a 3 away from the transition. 

C. Subleading IR divergences 

The analysis of Ref. 17] points out that, in addition to the leading order IR divergences, there are subleading 
IR divergent contributions. Although the former have been understood and summed using tm^PT, what of the 
latter? I have not done a complete analysis of this question, but I think the above analysis of the leading divergences 
makes the following conjecture reasonable: the subleading IR divergences are removed by a small shift in the vacuum, 
perturbatively calculable in tm^PT (except near the end-points). This incorporates the NLO and higher order effects 
of £5 and £q, as well as the contributions of higher order operators. In other words, I conjecture that there is no 
breakdown of the Symanzik expansion, in the sense described above, as long as one uses tm^PT to include the 
dominant terms (£5 in the GSM regime and £q in the Aoki regime) non-perturbatively. Progressively higher order 
terms (£7, £g, etc.) will then have progressively smaller effects, each subsequent order suppressed by oAqcd- 

This brings me to the logical relation between the approach of Ref. 01 an d tm^PT. I want to reiterate that the 
starting point of tm^PT is the Symanzik expansion, but that tm^PT goes beyond by including our non-perturbative 
understanding of spontaneous chiral symmetry breaking. Thus it is based on a tower of two local effective theories: 
the first describing the long-distance physics of lattice QCD, and the second describing the long-distance behavior of 
the first. The theoretical status of the second effective theory, i.e. chiral perturbation theory, is, perhaps, less solid 
than the first, since one cannot do a perturbative analysis following Symanzik as chiral symmetry breaking is a 
non-perturbative phenomenon. Nevertheless, the theoretical basis for tm^PT, which is identical to that for continuum 
%PT, is strong 19|,|20|. The main shortcoming of tm^PT, I believe, is not its theoretical foundation, but rather the 
practical need to truncate the expansion at NLO or, perhaps, NNLO, and the related issue of how well the expansion 
converges for present quark masses and lattice spacings. In this regard, the present analysis and the conjecture made 
above are important, because they imply that there is no non-uniformity in a joint expansion in m q and a. The 
absence of IR divergent terms implies that one makes only small errors by truncating the Symanzik expansion at, say, 
£7 (as is done in the previous section). 



19 



VI. CONCLUSIONS 

Studies of tmLQCD are still at a relatively early stage, and it is important to fully understand the practical and 
theoretical issues that are involved. Key issues are the determination of the critical mass and the size, and chiral 
behavior, of the remaining discretization errors. Particularly important discretization errors are those giving rise to 
parity and flavor breaking, for unless these are small in practice they will likely lead to complications in extracting 
phcnomcnologically important hadronic quantities. 

Here I have focused on the theoretical issues, showing how by applying chiral perturbation theory to Symanzik's 
effective Lagrangian one can systematically study the properties of discretization errors. The aim here is to understand 
the errors so that they can be removed, or reduced, in a controlled way. For example, if one wants to work at maximal 
twist, where the errors of 0(a) are automatically removed, then it is useful to know what happens when the choice of 
critical mass has errors of a given size. The forms given in Ref. |13| and generalized here provide such information to 
NLO in the joint chiral-continuum expansion. In particular, I have extended previous results in the Aoki regime to 
NLO, allowing consistent NLO fits for all quark masses. 

On a practical level, there is general agreement in the literature on how to proceed: one should use one of the 
"optimal" non-perturbative definitions of maximal twist [methods (i), (ii) or (iii)]. Here I have clarified the relation 
between the different definitions, both in the GSM and Aoki regimes. 

There is however, disagreement on the interpretation of what happens to discretization errors as one approaches 
the chiral limit. With non-optimal choices of the critical mass (including the a m 7T = 0" definition irrespective of the 
accuracy of the extrapolation) the results for a number of physical quantities exhibit a "bending" effect, in which 
the results diverge away from the expected chiral behavior for m q < a. These can be understood qualitatively in 
tmxPT as a result of working at a non-maximal twist angle, with the divergence from tt/2 increasing toward the chiral 
limit 13]. Examples of the expected forms are given in Figs. [S] and [fj] A seemingly different interpretation has been 
given in Ref. [T^ , in which the bending is due to discretization errors which are proportional to inverse powers of the 
quark mass, and which, for m q ~ a, arise from contributions proportional to all powers of a. I have shown, however, 
that the interpretations are consistent: the apparently IR divergent errors of Ref. lfj are simply an expansion of the 
convergent expressions of tnr^PT, and indicate the need to do a non-perturbative shift in the vacuum. 

A similar discussion applies to the discretization errors that arise when one works with an optimal choice of 
maximal twist. The IR divergent errors analyzed in Ref. [TtI ] are summed up by tm^PT applied in the Aoki regime. 
One consequence is that there is no barrier to working at maximal twist in this regime, i.e. with m q ~ a 2 , as long 
as one uses the optimal choice of maximal twist determined in this regime (rather than by extrapolation from larger 
masses) and as long as one does not work all the way down to the second-order end-points if one is in the scenario 
with a first-order transition. 

It is suggested in Ref. 0] that the apparent IR divergences signal a breakdown of the Symanzik expansion at small 
quark masses (m q ~ a 2 for an optimal choice of maximal twist). I have argued that this is not the case, although 
the situation is somewhat subtle. The apparent IR divergences are an indication that, in general, the leading terms 
in the Symanzik expansion (those proportional to a and a 2 ) must be treated non-perturbatively, and this is what 
is accomplished by tm%PT. This might be interpreted as a breakdown in the simple expansion in powers of lattice 
spacing, since insertions of any number of factors of i§ and £q are required. On the other hand, one can treat 
subsequent terms, £7, is, etc., as small corrections, so in this sense the expansion remains intact. In fact, at maximal 
twist, defined optimally, the contributions of £5 and £q can also be treated perturbatively, so the Symanzik expansion 
remains valid in the usual sense that all powers of a are manifest. It is for this reason that automatic improvement at 
maximal twist remains valid into the Aoki region (and all the way down to m q = in the Aoki-phase scenario [T^.ll3|). 
This result should simplify numerical applications of tmLQCD. 

After this paper was completed, Ref. [3(J appeared, containing results overlapping those in sees. ITT1 ITTT1 and llVl as 
well as detailed fits of numerical results to tm%PT expressions. 
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